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A study of hadronic data up to TEVATRON energies in terms of relativistic hydrody- 
namics indicates an extended 1-dimensional stage of the expansion which suggests a jet 
like behaviour of the fireball along the collision axis. 

The Landau Hydrodynamical Model (LHM) exists already for 40 years but there has 
been little progress in understanding its successes. Hydrodynamics assumes in general 
local thermal equilibrium (I.e.), a condition difficult to realize in small and short lived 
hadronic systems with a typical size of 10 -13 cm and a corresponding lifetime of ~ 10 -23 
sec. The discovery of subentities of hadrons (quarks and gluons) with the associated 
proliferation of degrees of freedom has facilitated the believe in I.e. M, nevertheless the 
phenomenological success of the LHM has not been understood so far. The situation is 
better for heavy ion reactions which are larger systems, and where it is easier to get I.e. 

We want to suggest that one possible reason for the success of LHM in hadronic reactions 
is the fact that Landau Q and most of his followers used only a 1-dimensional (Id) solution 
which corrects for the possible absence of I.e. in these reactions. This conclusion follows 
from a comparison of Id and 3c? solutions to be reported below. The Id approach assumes 
that the strongly compressed initial fireball expands at first mainly in the longitudinal 
direction (the width of the rapidity distribution is directly connected with the strength 
of this flow). Then a conical (3d) correction follows. A reasonable estimation for the 
moment r^a when the conical expansion starts (R denotes the transverse radius of the 
fireball) is: 



t-u = Vt 2 - x 2 = a 3d R (1) 
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where a^a is a phenomenological parameter (determined by Landau from simple geome- 
trical considerations [Q] to be equal to a 3 d = = (1 + Cq)/cq; c is the speed of sound). 

In Fig. 1 we present rapidity distributions calculated for "conical" Id and 3c? || solutions 
compared with pp data at SPS (y/s = 20 GeV) and ISR (y/s = 53 GeV) energies. The 3d 
solution uses Co as given by lattice QCD || and K =0.35 at y/s = 20 GeV and K = 0.176 
at y/s = 53 GeV. More realistic values of K lead even to a worsening of the agreement 
with data. The Id solution uses c\ = 0.18 at y/s = 20 GeV and Cq = 0.195 at y/s = 53 
GeV and K = 0.5. 

It turns out that only the Id solution is able to fit the data with reasonable values for 
Co and inelasticity K. The transverse expansion (present only in the 3d case) develops at 
the expense of the longitudinal expansion and therefore reduces the width of the rapidity 
distribution. However, going to still higher energies we have found that one has to increase 
the duration of the Id stage even further (by increasing a 3 d above the limit given by 
Landau). In Fig. 2 we show fits to different pseudorapidity distributions which are 
obtained with a^d increasing from y/2a L at y/s~ = 53 GeV to Vl0a L at 1800 GeV. We 
have checked that these results hold (almost) independently of the concrete variant of 
initial conditions and equation of state (EOS) provided they are physically reasonable. 

This observation poses a serious problem for LHM because the corresponding extended 
Id stage is not present in the "real" 3c? dynamics (which is based on the assumption 
that each fluid cell has in its rest frame an isotropic pressure - a result of the assumed 
isotropic momentum distribution corresponding to local equilibrium (I.e.)). One might 
therefore argue that because of the breakdown of I.e. conventional hydrodynamics is not 
valid anymore. On the other hand, the success of Id LHM in describing data, illustrated 
above, allows also a different interpretation. In the following we shall argue that at high 
energies a "new" physical effect occurs, namely a strong anisotropy in the flow caused by 
some physical processes acting on top of the conventional hydrodynamical description. 

The simplest "model" for such an anisotropy would be to postulate the existence of 
beam jets associated e.g. with the leading particles. This would make necessary a re- 
formulation of the inelasticity effect in the LHM. So far inelasticity K was considered in 
the LHM by assuming that only the function K of the available energy contributed to 
the mass of the fireball which underwent hydrodynamical expansion In this way 

the leading particles "had done their job" and did not interfere anymore with the cen- 
tral fireball. This treatment may be an oversimplification. A more realistic approach in 
this direction is represented by the two component model proposed in ||. It is based 
on the analysis of multiplicity distributions P(n) at energies between 20 and 540 GeV 
||. They were interpreted as indicating the presence of two different types of sources 
emitting secondaries: (i) - chaotic, provided by gluonic interaction (i.e., equilibrated) and 
concentrated in the central rapidity (i.e., hydrodynamical) region with P(n) of negative 
binomial type and (ii) - coherent, provided by the leading valence quarks (therefore far 
from equilibrium) and extending over the entire rapidity region (but contributing mainly 
to fragmentation region) with P(n) consistent with a Poisson distribution. In this context 
the anisotropy and the elongation of r 3 d can be viewed as a manifestation of the coherent 
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component due to the leading valence quarks|]. 

Formally these possibilities could be formulated in terms of anisotropic hydrodynamics 
as proposed in J7J. In the present study we shall limit ourselves just to consider the 
extended Id stage as a phenomenological observation which may have consequences for 
the interpretation of data from future experiments for hadronic and heavy ion collisions 
(RHIC, LHC or SSC). Here the hadronic reactions provide a lower limit for stopping, 
lifetime and equilibration and an upper limit for the width of the rapidity distribution a. 
For this last quantity we get as an upper limit 

a < ^fhiR/Si + atherm] (cr therm < 1.Q5 for 7/ < 0.2 GeV). (2) 

where <5j is the initial longitudinal extension of the fireball and o t h erm the contribution of 
the thermal motion to the rapidity width. It grows with energy slower than the phase 
space. This is shown in Fig. 3. A hydrodynamic stage in high energy collisions leads 
therefore to a limited value of a (e.g. a < 5 at LHC energies). It can be also shown that all 
relevant information concerning the transition from a strongly interacting non-equilibrium 
system to a thermalized fireball is contained in the fragmentation region (i.e., the phase 
space region where the transition from a local equilibrium in a pre-equilibrium stage takes 
place. This is also illustrated in Fig. 4 where one can see how the initial longitudinal size 
of the fireball (determined in the pre-equilibrium stage) is strongly reflected in the shape 
of rapidity distribution. The main effects appear in the fragmentation region. 

We conclude than that a hydrodynamical analysis of pp data indicates a large extension 
of the Id stage of the expansion and is described approximately by the Khalatnikov 
solution [0. The observation of the fragmentation region (3 < y cm < 5) is essential for 
the investigation of the transition from the pre-equlibrium to the local equilibrium stage 
of the reaction. 
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